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ABSTRACT 

Context. With their period-luminosity relation, "Classical Cepheids" (CC) are the most common primary distance indicators within 
the Local Group, also providing an absolute calibration of important secondary distance indicators. However, the predicted position 
of these pulsators in the HR diagram, along the so called blue loop, that is the expected distribution of Cepheids within the instability 
strip is affected by several model inputs, reflecting upon the predicted PL relation. 

Aims. The aim of this work is to quantitatively evaluate the effects on the theoretical PL relation of current uncertainties on the chem- 
ical abundances of Cepheids in the Large Magellanic Cloud (LMC) and on several physical assumptions adopted in the evolutionary 
models. We will separately analyse how the different factors influence the evolutionary and pulsational observables and the resulting 
PL relation. 

Methods. To achieve this goal we computed new sets of updated evolutionary and pulsational models. 

Results. As a result, we find that present uncertainties on the most relevant H and He burning reaction rates do not influence in a 
relevant way the loop extension in temperature. On the contrary, current uncertainties on the LMC chemical composition significantly 
affect the loop extension and also reflect in the morphology of the instability strip; however their influence on the predicted pulsa- 
tional parameters is negligible. We also discussed how overshooting and mass loss, sometimes suggested as possible solutions for the 
long-standing problem of the Cepheid mass discrepancy, influence the ML relation and the pulsational parameters. 
Conclusions. In summary, the present uncertainties on the physical inputs adopted in the evolutionary codes and in the LMC chemical 
composition are negligible for the prediction of the main pulsational properties. On the other hand, the inclusion of overshooting in 
the previous hydrogen burning phase and/or of mass loss is expected to significantly change the resulting theoretical pulsational sce- 
nario for Cepheids, as well as the calibration of their distance scale. These systematic effects are expected to influence the theoretical 
Cepheid calibration of the secondary distance indicators and in turn the resulting evaluation of the Hubble constant. 

Key words, cepheids - stars:evolution - stars: distances 
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1. Introduction 

During the central helium burning phase, intermediate mass stars 
show an excursion toward higher effective temperature, in the 
Hertzprung-Russell (HR) diagram, subsequently coming back 
toward the asymptotic giant branch (blue loop). During this 
phase the stars cross the instability strip, becoming Classical 
Cepheids. These variable stars, thanks to their characteristic 
Period-Luminosity (PL) relation, are the most common pri- 
mary distance indicators within the Local Group. Moreover, 
from space, they are observable al so in external galaxie s, up 
to distances of 20-30 Mpc (see e.g. [Freedman et al. 2001). On 
this basis they provide the absolute calibration of important 
secondary distance indicators, such as the maximum luminos- 
ity of Supernovae la, the Tully-Fisher relation, surface bright- 
ness fluctuations, and the planetary nebulae luminosity func- 
tion, that are the basis for measurement of the Hubble con- 
stant (Freedm an et alJl200U ISaha et alj|2001l) . Any systematic 
error affecting the Cepheid PL relations affects the extragalac- 
tic distance scale and the estimate of the Hubble constant. 
Therefore, even though the Cepheid PL relation is assumed to 
be universal and calibrated throu gh the sample in the Large 
Magellanic Cloud (LMC, see e.g. iMadore & Freedmanl fl991; 



lUdalski et al Jl 19 993, the possibility that the Cepheid properties 
and PL relation depend on the chemical composition of the par- 



(see e.g. Sakai et alj 2004 


; Macri et alj|2006l iRomaniello et al.l 


2005, 2008; Marconi et al. 


l2005tlBono et all2008l). On the other 



hand the theoretical prediction for the existence of a PL rela- 
tion for Classical Cepheids relies on the assumption that inter- 
mediate mass stars undergoing central helium burning are char- 
acterized by a mass-luminosity (ML) relation, as predicted by 
stellar evolution models. On very general grounds the pulsa- 
tion period is expected to be tightly correlated with the lumi- 
nosity, the mass and the effective temperature of the star, for 
each given chemical composition. It is only thanks to the ML 
relation that the period can be simply related to the luminosity 
and the effective temperature, allowing a theoretic al prediction 
for the period-luminosity-color relation (PLC, s ee Feast 1984 
199 11 [1995 ICaldwell & Coulsod|198d lLanev & Stobiejll98e 
lstotherslll988t ICaldwell & LanevT " 



1991b ISasselov et al. 



19971) . 
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However the theoretical ML relation depends on the several 
physical processes included in stellar evolutionary codes as well 
as on the adopted chemical composition. 

By averaging the PLC relation on the color extension of 
the pulsation instability strip one obtains the PL relation. This 
implies that the PL relation has a statistical nature and ob- 
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servationally holds only when a statistically significant sample 
of Cepheids (at the same distance) is available. Thus the pre- 
dicted temperature extension of the blue loop for a given mass 
also reflects upon the expected distribution of Cepheids within 
the instability strip and in turn upon the theoretical PL rela- 
tion. This is a very critical point because the predicted exten- 
sion of the blue loop in the HR diagram is affected by several 
model inputs (chemical composition, efficiency of external con- 
vection, nuclear burning cross sections, etc..) producing effects 
that are not yet completely understood (see e.g. lXu & Lil l2004a; 
IStothers & Chinll994tlBrunish et al.ll990t) . In this context, a de- 
tailed comparison of evolutionary and pulsational models with 
the observations can provide insight into the physics of stellar 
evolution and pulsation. 

The aim of this work is to quantitatively evaluate the ef- 
fects on the predicted PL relation of current uncertainties on 
the chemical abundances of Cepheids in the LMC and of sev- 
eral physical assumptions adopted in the evolutionary models. 
In particular, we will separately analyse how the different factors 
influence the evolutionary and pulsational observables (exten- 
sion in effective temperature of the blue loop, Mass-Luminosity 
relation, morphology of the instability strip, period distribution) 
and the resulting PL relation. For the chemical composition we 
choose that of the LMC with its associated uncertainty. All the 
relations based on the newly computed sets of evolutionary and 
pulsation models will be made available to the community. In 
a forthcoming paper (Valle et al. 2009 in preparation) we will 
focus on the comparison of these new predictions with the ob- 
servations. The organization of the paper is the following: in 
Sect. |2] we present the theoretical scenario both from the evo- 
lutionary and the pulsational point of view; Sect. [3] contains a 
general discussion about the problems in calculating blue loop 
models; in Sect[4]the dependence of the predicted extension of 
the blue loop on the uncertainties affecting nuclear reaction rates 
is investigated; in Sect. [5] we focus on the effects of chemical 
abundance uncertainties on the most relevant evolutionary and 
pulsational features, whereas in Sect. 16.11 we take into account 
possible noncanonical effects in the evolutionary scenario. The 
conclusions, reported in Sect. [7] close the paper. 

2. Theoretical scenario 

In this section we present the new sets of evolutionary and pul- 
sation models that have been used for the calculations. 



2.1. Evolutionary stellar models 

Evolutionary tracks were computed with an updated version of 
the FRANEC evolutionary code (see e.g. Chieffi & Str anierol 
Il989t iDeglTnnocenti et all l2008h including the O PAL 2006 
equation of state (EOS)Q (see also Rogers et al. 1996) and radia- 
tive opacity tabled (see also Iglesias & Rogers 1996) for tem- 
peratures higher than 12000 K. Moreover our models i nclude the 
conductive opa cities by IShternin & Yakovlevl d2006l) (see also 
Potekhin 1999) and the atmospheric opacities bv lFerguson et alj 
(l2005l) PI. All the adopted opacity tables have bee n calculated by 
assum ing the solar relative metal abundances by lAsplund et alj 
(120051) . 

Since the 12 C(avy) 16 and 14 N(p,y) ls O reaction rates af- 
fect the evolutionary properties of intermediate-mass stars 



and in particular the t emperature blue loop extension (see 
e.g. iBrunish et al] Il990h . we adopted the updated rates by 
Ham mer et alj d2005l) for the l2 C+ a and by the LUNA col- 
laboration (see Imbriani et al. 2005, and references therein) for 
the 14 N+p cross sections. F or the other nuclea r reactions we 
adopted the same rates as in Ciacio et al. d 1 997b . The code im- 
plements weak screening (Salpeter 1954) , weak-in termediate, 
intermedi ate- strong screening (jGraboske et al.l 1973b . and strong 
screening dltoh et alJI 197711 1979b . 

For convective mixing, we adopt the Schwarzschild crite- 
rion to define regions in whic h convection elemen ts are accel- 
erated (see the description in iBrocato et alj|2003l) with an ad- 
ditional algorithm to take into account the in duced semicon- 
vection during the central He-burning phase (Castellani et alJ 
1 1971 alibi). Rotationally ind uced effects (see e.g. Maeder & Zahn 
1998; Pala cios et al.ll2003l) are not included in our models. 

To model envelope co nvection we adopted, as usual, the 
mixing length formalism (Boh m^Vitensd fl958l) . The mixing 
length parameter, a, governing the efficiency of convection, has 
been calibrated to reproduce the observed stellar colors in the 
Large Magellanic C loud cluster NGC 1866, obtaining a=1.9 
dBrocato et alJ l2004). This result is also dependent on the atmo- 
spheric models adopted to transform evolutionary calculations 
from the theoretical (logL-logT c ) to the observational plane. In 
this paper w e adop ted the color-temperature transformations by 
lBesseletalJdl998l) . 

The original helium mass fraction (Y) cannot be measured 
spectroscopically and it is generally estimated by assuming a 
linear relation between Y and Z: 
AY 

Y = Y P + — x Z (1) 
AZ W 

where Yp is the primordial helium abundance and AT/AZ is the 
estimated helium-to-met al enrichment ratio. For Yp we adopted 
the value Y P = 0.248 dPeimbert et alj|2007l; llzotov et alj|2007h 
and AT/AZ = 2 dJimenez et alj|2003l) . 

The metallicity adopted for the models is consistent with 
the most recent estimates of [Fe/H] for LMC Cepheids (see 
iLuck et alj|1998t iRomaniello et alj|2008h and also with several 
independent determinations of the LM C metallicity (see e.g. 
iKorn et alj|2000l lAndrievskv et alj|2001l) . 

The total metallicity can be expressed as follows: 



AY 1 

Z = (1 - Yp) 1 + — + 

' ' AZ (ZX) e 



\Q-[FelH] 



(2) 



1 http://www-phys.llnl.gov/Research/OPAL/EOS.2005 / 

2 http://www-phys.llnl.gov/Research/OPAL/opal.html 
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available at http://webs.wichita.edu/physics/opacity 



where (Z/X) Q is the solar metal to hy drogen abundance rat io. 

In this work we adopt the updated Aspl und et al.l {2005 ) mix- 
ture, based on the analysis of spectroscopic data by means of 
tridim ensional hydrodynamical atmospheric models, instead of 
that bv lGrevesse & Noelsl (119931) . hereinafter GN93, widely used 
in the past. This leads to a significant reduction of the the Z/X 
ratio for the Sun, from Z/X * 0.0245 to Z/X « 0.0165. As a 
result the traditionally assumed metallicity for LMC decreases 
from Z « 0.008 to Z « 0.005. With this choice of metallicity, 
from Eq.Q] one obtains Y ~ 0.26, which we assume as our stan- 
dard value. 

Several investigations have dealt with the analysis of the ef- 
fects of this revision of the solar abundance, revealing discrepan- 
cies between p redicted and helioseismologically infer red quan- 
tities (see e.g. lBahcall et al.ll2005l;lBasu & Antiall2004l) . This sit- 
uation is further complicated by the fact that current uncertain- 
ties on the physical inputs adopted in the solar models can af- 
fect the results of the quoted investigations. Moreover some au- 
thors questioned the accuracy of the oxygen abundance derived 
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bv lAsplund et al.l(l2005l) (see e g.lCaffau et al. 2008). For these 
reasons, even if we adopted the Asplund et al. (2005) solar mix- 
ture as our reference mixture, in the following we will take into 
account the uncertainty related to this assumption and its influ- 
ence on the LMC metallicity. 

2.2. Nonlinear pulsation models 

Our pulsational models are based on a nonlinear non- 
local time-dependent convective hy d rodyn amical code 
dBono & Stellingweri Il994t iBono et al l Il999al and refer- 
ences therein). This code has been suc cessfully appl ied to 



taD 



2003 



2004 



several classes of pu l sating stars (see e.g . [Marconi et 
Pi Criscienzol 120041: IBono et al] tl 999b; Marconi et 
2005: [Marconi & Di Criscienzol 120071 and references therein). 
In this paper the physical ingredients of Cepheid pulsation 
models, and in particular the stellar opacity tables, have been 
updated in order to insure consistency between the pulsation and 
the evolutionary models. The new pulsation models have been 
constructed with the same chemical abundances of evolutionary 
models, covering the same range in stellar mass and assuming 
the luminosity p redicted by the c anonical Mass-Luminosity 
(ML) relation bv lBono et ail d2000l) . For Z = 0.005, a second 
luminosity brighter by 0.25 dex has been included to account for 
the effect of mild overshooting and mass loss (see Subsec. l6Tl ). 



3. Modeling the blue loop morphology 

The extension in temperature of the blue loop is an impor- 
tant evolutionary feature. Stellar models show that the excur- 
sion to the blue is p roduced by an envelope contraction (see e.g. 
iRenzini et alJI 19921) : models with more extended loops undergo 
a more significant envelope contraction which heats the hydro- 
gen shell, increasing its luminosity. The ret urn of the track to th e 
red is due to a variety of effects (see e.g. ISestito et al.l 12002) : 
a) when the central helium abundance decreases below Y c ~ 0.1 
the star reacts to the decreased burning rate by means of a core 
contraction increasing again the energy production and thus the 
amount of energy reaching the envelope, which expands and 
cools; b) the H shell reaches the discontinuity in the chemical 
profile left by the first dredge up before the core contraction (see 
e.g. the discussion in lLauterborn et al.l jl97lt iHuang & Weigert 
1983; ISestito et alj|2002t IXu & Lill2004al). However, as widely 
discussed in the literature (see e.g. IStot hers & Chin 1993; 
Alcock & Paczvnsldl978l:lRobertsonll97l[:lBrunish et al.ll 990; 



Sestit o et al.ll2002h IXu & Lill2004bl) . the formation and the ex- 



tension of loops appear to be caused by the interaction of several 
factors and is extremely sensitive to small changes either in the 
stellar chemical composition or in the physical inputs adopted in 
the calculations. Several authors have discussed the poor stabil- 
ity naturally existing in many stellar models which lie on the red 
supergiant branch during the central He burnin g phase, that is, 
around the beginning of the b lue loop (see e.g. Noels & Gabriel 
1 19731: IStothers i & Chinlll991al) . 

The dependence of the loop morphology on changes in the 
chemical composition will be analysed in Sec. [5] In general, an 
increase of the original He abundance increases the mean molec- 
ular weight while it leads to a decrease of the mean opacity (see 
e.g. IVemurv & Stotherslll978l) . even if the dependence on He 
abundance is not linear with metallicity. On the other hand a 
metallicity variation changes the blue loop extension ma inly due 
to the opacity variation (see e.g. Stot hers & Chin|[l993l) even if 
the behaviour can be more complicated and a quantitative expla- 



nation of the loop dependence on chemical composition is still 
lacking in the literature. 

The evolutionary properties of intermediate and high mass 
stars in the central He burning phase strongly depen d on the ef- 
ficiency of the He burning nuclear rea ctions (see e.g. Ibenl l 19721: 
Brun ish et alll99CtlBono et alj2000h . This point is discussed in 
the following subsections in the light of updated nuclear cross 
sections measurements. 

Regarding the opacity dependence, any increase of the opaci- 
ties deep in the stellar interiors could lead to fainter luminosities, 
cooler effective temperatures and longer lifetimes. Regarding 
the loop extension, in principle, a higher opacity leads to a re- 
duc ed blue loop but the e ffect is in general very small (see 
e.g. IStothers & Chinlll991bl) . In addition, updated opacity tables 
agree amongst each other within 5-10% (see e.g. dis cussion in 
iBadnell et al.ll2005t iNeuforge-Verheecke et al.l 12001) thus very 
small opacity uncertainties are expected. Moreover, molecular 
opacities hav e almost negligible eff ects on the blue loop exten- 
sion (see e.g. lStothers&Crrinlll993l). 

Several authors (see e.g.lChiosi et ailll992t IStothers & Chinl 
1992: iBrocato & Castellani| [l993) noticed that the extension of 
the convective H burning core affects the following He burning 
phase. This involves the convective core overshooting efficiency, 
whose influence on the blue loop extension will be discussed in 
detail in Sec. 16. II and the method adopted for treating the semi- 
convective layers that develop around the H burning core close 
to the central H exhaustion. If a significant amount of overshoot- 
ing is included in the calculations, during the core He burning 
phase the luminosity is brighter, the lifetime is shorter and the 
blue loop is less extended and rises less steeply than in absence 
of overshooting (see e .g. iMatraka. Wassermann. Weigertl [l982; 
IChin & Stothers|[l99Th . 

The possible occurrence of convective shells and the 
consequent modification of the internal chemical profile 
are strictly related to the adopt ed stability criterion, either 
Schwarzschild or Ledoux (see e.g.[Salasnich. Bressan & Chiosil 
1999; Stothers & Chir3ll994tlCanutoF oOO). However the stabil- 
ity criterion has a marginal effect on stellar models ranging from 
»3M to «T2 Mp, as the ones taken into account in our work 
IStothers & Chinl d 19931) . 

Moreover downward convecting overshooting from the outer 
convection zone which modifies the envelope chemical profile 
influences the formation of blue loops. 

For the same reason the occurrence of blue loops is thus in- 
fluenced by the mixing length value a, which governs the con- 
vection efficiency in the external convective envelope. However 
when a blue loop develops (it is not true in some cases of small 
and moderate a values) the maximum effective temperature dur- 
ing the blue l oop was found to be o nly a very weak function 
of a (see e.g. IChin & Stofhersl[l99ll) . Moreover, since the at- 
mospheric model is fixed, a is generally calibrated through the 
comparison between theory and observation for the stellar color. 
For this reason, in this paper we will not include variations of 
this parameter, which is fixed to a=1.9. 

Another important phenomenon affecting the evolution of 
He burnin g intermediate mass star s is mass loss (see e.g. the 
rev iew by |Chi osi & Ma ederl 1 19861) or the more recent work 
by [Salasnich, Bress an & Chiosil d 1999b. As discussed b y sev- 
eral authors (see e.g. Salasnich, Bress an & Chiosil [T999) mass 
loss mainly influences the evolution of masses higher than a; 15 
M Q but its influence on the loop morphology, as discussed in 
Sec. 16.21 can be significant also for lower masses. However 
the physical mechanisms which affect mass loss and its depen- 
dence on metallicity are not still completely understood and 
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thus the range of possible efficiencie s for the phenomenon is 
very wide (see e.g. the discussion in Mae der & Mevnell F2000; 
Vazq uez et al.l 120051) . In this paper we work in the canonical 
regime, also because we concentrate on intermediate mass stars 
but to take into account a possible mass loss contribution, in 
Sec. 16. 21 we discuss some additional models computed with dif- 
ferent mass loss efficiencies. 

Due to the dependence of the loop extension on several phys- 
ical parameters, if one compares recent sets of intermediate mass 
stars models available in the literature (which in general differ 
by the assumed chemical composition, the adopted physical in- 
puts, the external convection efficiency etc..) one finds signifi- 
cant variations in the loop morphology. In the following sections 
we will discuss the main parameters which have already been 
demonstrated to potentially affect in a relevant way the blue loop 
morphology: He burning cross sections, chemical composition, 
overshooting and mass loss efficiency. 



4. Dependence of the loop extension on the 
uncertainties on the nuclear reaction cross 
sections 

Several authors have investigated the dependence of the loop ex- 
tension on the H burning (thr ough the CNO cycle) and the He 
burning reaction rates (see e.g.lBrunish et alJl990HRenzini et all 
ll992tlSesfito et al.l2002l:IXu & Lil2004atlWeiss et alj2005l) . and 
have shown that the loop morphology is sensitive to variations 
of the adopted cross sections. 

In this section we focus on the dependence of the loop forma- 
tion and its temperature extension on variations of the relevant 
nuclear fusion cross sections within their estimated uncertain- 
ties. For stars in the central He burning phase, the most relevant 
nuclear reactions are the 3a, the 12 C(a,y) 15 0, that affects the He 
burning efficiency and the 14 N(p,y) 15 that governs the shell H 
burning through the CNO cycle. 

4.1. l4 N(p,y) l5 cross section 

In stars with internal temperatures higher than ss 2 • 10 7 K, the 
hydrogen burning is dominated by the CNO cycle whose rate 
is determined by the slowest process, the 14 N(p,y) ls O reaction. 
For this reaction the Gamow peak is at « 30 keV; thus the cap- 
ture cross section must be known down to this energy at which 
it has a very low value, preventing a direct laboratory measure- 
ment. Instead, the cross sections are measured at higher ener- 
gies and then extrapolated, by means of the astrophysical S(E) 
factor, down to the Gamow peak region. The occurrence of res- 
onance and resonance tails complicates the extrapolation and it 
results in very large uncertainties on the evaluated cross sections. 
Recently, the 14 N(p,y) I5 cross section has bee n measured by 
the LUNA collaboration close to 70 keV (see e.g . llmbriani et all 
l2005t iBemmereret alJl2006t iLemut et al.1l2006h . leading to an 
improved estimate of the S(E) factor at astrophysical energies 
(5(0) = 1.61 ± 0.08 keV, about half of the value of the NACRE 
compilation by lAngulo et alJ (Il999h t. 

To quantify the effect of a variation of the I4 N+p cross sec- 
tion, we computed stellar models for several stellar masses, vary- 
ing the 14 N+p reaction rate of the estimated uncertainty with 
respect to the LUNA value (adopted in our standard models), 
models IX and X of Table Q] For all the investigated cases the 
effect on the loop morphology of the present uncertainty on the 
14 N(p,y) 15 rate is negligible, suggesting that this uncertainty is 



not expected to contribute significantly to the final error budget 
for the predicted Cepheid pulsation properties. 

4.2. 3a cross section 

Recently Fvn bo et ail (12006) reported new measurements for the 
3a rate, finding, in particular, a resonance at energies of «11 
MeV; however, in the temperature range of interest, the differ- 
ences with respect to the most widely adopt ed 3a rate are within 
the uncertainty evaluated by the NACRE (Ang ulo et all fl 999) 
compilation, reachin g a maximum of ^ 20% at temperatures of 
about 10 8 K (see e.g. lWeiss et alJl2005l) . 

For this reason we decided to take the error quoted by 
NACRE (20%), in the temperature range of interest, as a con- 
servative estimate of the uncertainty on the 3a reaction rate. We 
computed specific models by varying the 3c rate within the esti- 
mated uncertainty for different stellar masses (in our "standard" 
case) finding no relevant variations in the loop extension. These 
results suggest that the present uncertainties on the 3a rate do not 
substantially influence the evolutionary and pulsational quanti- 
ties for Cepheids stars. 

4.3. 12 C(a,y) 16 cross section 

The measurement of the 12 C + a cross section at the energies of 
interest in stars 300 keV) is presently impossible due to the 
smallness of the predicted value (=; 10~ 17 b). The extrapolation 
from higher energies is quite difficult, because the cross section 
in this energy region is a mixture of groundstate and cascade 
transitions, and the cross section for the ground state transitions 
is dominated by the tails of subthresh old resonances with th e 
interference of other processes (see e.g. Buchm ann et al.l l2006). 
For these reasons the uncertainty quoted in the past for this rate 
was very large, up to a factor of two, and the ef fect of this er- 
ror on the loop extension was very strong (see e.g. Brunish et al. 
Il990l) . Given the astrophysical importance of this reaction, in the 
past decade a significant experimental effort has been devoted to 
the impro vement of the determ ination of its cross section. Recent 
results by [Hammer et al. | d2005|) of direct 12 C(a,y) 16 measure- 
ments in the temperature range 0.001 < Tg < 10, Tc, being the 
temperature in 10 9 K, lead to a new determination of the extrap- 
olated astrophysical reaction rate with a maximum total uncer- 
tainty of + 25% (see lHammer et al.ll2005l) . 

In the present work we adopt for the 12 C(ar,y) 16 reac- 
tio n rate the parametriz ation and the parameter values given 
bv lHammer et al.1 d2005): at the typical He burning temperature 
(10 8 -h2-10 8 K) this new rate is about 20% lower than the NACRE 
rate. 

As discussed in iBrunish et al.l (119901) . the blue loop exten- 
sion depends on the I2 C + a cross section: the higher the rate of 
this reaction, the higher the effecti ve temperature excursion of 
the loo p. Several authors (see e.g. ISestito et alJ 120021: IXu & Lil 
2004a) showed that the 12 C(a,y) 16 cross section is related to 
the efficiency of the envelope contraction and thus to the loop 
extension. 

We evaluated the influence on the loop extension of the present 
uncertainty on the l2 C{a,y) l6 reaction rate for models of dif- 
ferent masses. To this aim we calculated two additional sets of 
tracks with "standard" c hemical compo s ition ( Z=0.005 Y=0.26) 
and the 12 C+a rate by lHammer et al.l (120051) increased or de- 
creased by the estimated uncertainty of 25% (models VII and 
VIII of Table[T|i. As shown in Fig.Q] available online, the effect 
is negligible for all the masses. 
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Thus, at variance with results published in the past, now the 
uncertainty in the 12 C+ a nuclear rate is noticeably reduced 
so that its effect on the loop morphology appears negligible. 
Consequently, in this paper, we neglect the uncertainty associ- 
ated with this cross section in the evaluation of the error budget. 

In conclusion, if one takes at their face value the present es- 
timated uncertainties of the cross sections relevant for Cepheids 
star evolution, one finds that now the precision of the experi- 
mental measurements, at the energies of astrophysical interest, 
is high enough to rule out significant effects of these uncertain- 
ties on the loop extension and thus on the pulsational quantities. 

5. Effects of the chemical composition 
uncertainties 

In this section we present the effects of the LMC Cepheid chem- 
ical composition uncertainties on evolutionary properties such 
as the He burning loop extension and the ML relation, and on 
pulsation properties such as the instability strip morphology and 
the PL relations. 

For the LMC iron content we adopt [Fe/H] =-0. 40 that, as 
quoted abo ve, is consistent with th e estimates by 
19981) and iRomaniello et alj d2008l) . In particular 
1998 ) suggest a range of -0.55<[Fe/H]<-0.19. By adopting 



Table 1. Computed s ets of evolutionary tracks and pulsa- 
tion models GS98 = IGrevesse & SauvaS (fl998h . AGS05 = 
Aspl und et all d2005l) . 



Luck et al 



Luck et al 



the lAsplund et al.l (|2005) solar mixture, that is a (Z/X) Q value 
of 0.0165, from the assumed [Fe/H] one obtains (see Eq. |2]i 
Z = 0.005 (our "sta ndard" value). On t he basis of the observed 
range suggested by iLuck et all dl998l) we allowed Z to vary 
in the range 0.0035< Z <0.008. Taki ng into account the cur- 
rent uncertainties o n AT/AZ (see e.g. Pagel & Por tinaril l l998t 
iGennaro et al. 2008), we allowed a variation of the original he- 
lium content from 7=0.26 to 7=0.28. With Z = 0.005, Y = 0.26 
corresponds to AT/AZ = 2, while Y = 0.28 implies AT/AZ = 6. 

One should also take into account the uncertainty on the 
solar mixture. As already discuss ed, tridimensional h ydro- 
dynamical atmospheric models by lAsplund et alj d2005l) re- 
duced the derived abundances of CNO an d other heavy ele- 
ments with respect to previous estimates dGrevesse & Sauval 
1998, hereafter GS98). GS 98 already improved the mixture 
by IGrevesse & Noelsl ( 1993b . widely adopted in the literature, 
mainly revising the CNO and Ne abundance and confirming the 
very good agreement between the new photospheric and mete- 
oric results for the solar iron abundance. 

The change of the adopted heavy element ratios might affect 
the evolutionary calculations of LMC Cepheids in two differ- 
ent ways: 1) changing the evolutionary characteristics at fixed 
metallicity; 2) changing the LMC inferred metallicity from the 
observed [Fe/H]. Regarding the second point, the adoption of 
the GS98 solar mixture leads to a LMC metallicity of Z »0.008, 
which is within the previous quoted metallicity uncertainty. The 
first point is analysed by taking into account the effects on the 
evolutionary features of a solar mixture variation at fixed metal- 
licity, as due to its effect on opacity, nuclear burning and equa- 
tion of state calculations. 

In the following we present different sets of models vary- 
ing the chemical composition within the quoted ranges of un- 
certainty. The combinations of Y and Z values adopted for the 
different sets of models are summarized in TableQ] The first col- 
umn of Table [JJ gives the reference name of the treated cases 
(std, Zi ow , Zhigh etc.), the second and the third columns give the 
adopted Y and Z values, whereas the fourth column reports the 
adopted chemical mixture at fixed Z. The cases XI and XII will 
be presented in detail in Subsec lSTTl in which we will discuss the 
effect of a mild core overshooting in the H-burning phase. 



Case 


Y 


Z 


Mixture 


I (std) 


0.26 


0.005 


AGS05 


II (Z,„) 


0.26 


0.0035 


AGS05 


HI (Z high ) 


0.26 


0.008 


AGS05 


IV (Y high ) 


0.28 


0.005 


AGS05 


V (Y hit ,i„ Zfogfi) 


0.28 


0.008 


AGS05 


VI (Z high , GS98) 


0.26 


0.008 


GS98 


VII( 12 C(a,7) lb O),„ s , I 


0.26 


0.005 


AGS05 


VIII ( l2 C(a,y) 16 0) tow 


0.26 


0.005 


AGS05 


ix ( 14 N( P ,r)' 5 o) fasA 


0.26 


0.005 


AGS05 


X ( 14 N(p, 7 ) 15 0), mv 


0.26 


0.005 


AGS05 


XI(/ OV = 0.1H P ) 


0.26 


0.005 


AGS05 


XII (Z ov = 0.25/f p ) 


0.28 


0.005 


AGS05 


XIII (/ ov = 0.25/fp) 


0.26 


0.005 


AGS05 



For each case of Table[TJwe computed evolutionary and non- 
linear convective pulsation models spanning a stellar mass range 
from 3-4M to 13-14M , which covers the range of masses 
crossing the instability strip (see for details Sec |5.3t . 
All the computed sets for the cases I to VI of Table[JJare plotted 
in Fig. |4l the fits of the corresponding instability strips calculated 
as explained in Subsec |3T2l are superimposed on the evolutionary 
models. 

Some models show a smaller secondary loop during t he redward 
excursion in the central He burning pha se (see e.g. Bon o et ail 
2000; Salasnich, Bressan & Chios iil999l too), which is not a nu- 
merical artifact. The physical reason for this secondary loop is a 
slight increase in the size of the convective core, due to the opac- 
ity increase in the C and O rich nucleus, leading to small varia- 
tions of the H and He burning efficiency. The mechanism is the 
analo g of the breathing pulses occurring in low mass stars (s ee 
e.g. lSweigart & Demarque1ll972i[T97llCastellani et al.lll985l) . 

5.1. Loop extension 

Figure|2]shows, for a 7 M model, the effects on the blue loop ex- 
tension of variations in the helium and metal contents within the 
uncertainties of these quantities (cases from I to IV of Table [TJ. 
The upper panel shows that increasing the metal content, at fixed 
helium abundance, the extension of the blue loop decreases, 
whereas the lower panel shows that an increase of the helium 
content within its uncertainty range, at fixed metallicity, leads to 
more extended loops, even if the effect is reduced compared to 
that of lowering Z. The behaviour is similar for the other anal- 
ysed model masses. 

A simultaneous variation of Z and Y within their uncertainty 
ranges does not result in an additive effect on the temperature 
extension of the loop. An example of this behaviour is shown in 
Fig- 0(f), for which the set of evolutionary tracks was computed 
with the upper values of Y and Z, namely T=0.28 and Z=0.008 
(case V of Table[TJ. Beside the apparent non additivity in log T e , 
there is rather complicated mass dependence with a more regular 
trend for the intermediate mass values (=; 5 - 8M Q ) and an in- 
verted trend for higher masses. Thus one should be very careful 
to predict the extension of the loop when more than one param- 
eter is varied, as it cannot be simply evaluated by combining the 
results obtained for the variation of each single parameter sepa- 
rately. 

The effect of changing the mixture at fixed Z is shown 
by comparing Fig. |4] (c) (case III of Table [TJ, in which the 
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Table 4. The same as Table 2 but for the first overtone bound- 
aries. 



Fig. 2. Upper panel: effect of a metallicity variation on the loop 
extension for a 7 M ; the three different lines correspond to cases 
I,n,m of Table [TJ std case (Z=0.005, solid line) Z tow =0.0035 
(dashed line) and Z/,,^/, =0.008 (dot-dashed line). Lower panel: as 
in the upper panel but for a helium abundance variation and for 
cases I (std, F=0.26, solid line) and IV (Yj,;„/,=0.28, dashed line) 
of Table 1. 



set of tracks is computed with Z = 0.008 and the mixture of 
AGS05, with Fig. g] (e) (case VI of Table [TJ, computed with 
the same total metallicity (Z = 0.008) but with the GS98 mix- 
ture. To analyse this case we changed the original element abun- 
dances of our models and we calculated suitable opacity ta- 
bles making advantage of the codes available at the OPAL web 
siteQ; the EOS is also varie d following the mixture change (see 
DeglTnno centi et alj|2006l for a more detailed analysis of the 
effects of a mixture update on the stellar models). 
The mixture change at fixed Z has a negligible effect on the loop 
exten sion for M < 8M . For higher masses we find that models 
with lAsplund et alj d2005l) com position do not show blue loops 
at all, whereas the ones with the lGrevesse & Sauvall (119981) mix- 
ture display extended blue loops. As this effect is evident only 
for about half of the analyzed mass range we decided to neglect 
it in the following. 



5.2. Instability strip 

For each chemical composition, mass and luminosity, the modal 
stability of pulsation models has been investigated for both the 
fundamental and the first overtone mode. The effective tem- 
peratures (hereafter in K) of the hottest and coolest pulsating 
fundamental models, increased and decreased, respectively, by 



Case 


a 






b 




a 


First overtone blue edge 


I (std) 


3.939 


+ 


0.019 


-0.043 


± 0.006 


0.005 


II (Zw) 


3.942 


+ 


0.019 


-0.044 


± 0.006 


0.005 


in (Z m ) 


3.913 


+ 


0.046 


-0.038 


± 0.015 


0.013 


IV (Y high ) 


3.922 


+ 


0.027 


-0.038 


± 0.008 


0.007 


V (Y high , Z high ) 


3.954 


+ 


0.012 


-0.050 


± 0.004 


0.003 


VI (Z Ugh , GS98) 


3.962 


+ 


0.012 


-0.052 


± 0.004 


0.005 


First overtone red edge 


I (std) 


3.847 


+ 


0.010 


-0.023 


± 0.003 


0.003 


II (Z lm ) 


3.834 


+ 


0.004 


-0.018 


± 0.001 


0.001 


HI (Z Mgh ) 


3.840 


+ 


0.009 


-0.024 


± 0.003 


0.003 


IV (Yhigh) 


3.834 


+ 


0.004 


-0.018 


± 0.001 


0.001 


V (Yugh, Zhi g h) 


3.826 


+ 


0.004 


-0.018 


± 0.001 


0.001 


VI (Z h , gh , GS98) 


3.831 


+ 


0.014 


-0.021 


± 0.005 


0.006 



http://www-phys.llnl.gov/Research/OPAL/opal.html 



the extension of the pass in temperature adopted in our anal- 
ysis, 50 K, correspond to the fundamental blue and red edges 
(FBE and FRE). Similarly the first overtone blue (FOBE) and 
red (FORE) edges are evaluated. Table [3] reports the quadratic 
fit coefficients of the FBE and FRE, plotted in Fig. [4] as solid 
lines. The linear fits of the FBE and FRE are available on line in 
Table [2] whereas the linear fits for the FOBE and FORE, plot- 
ted in Fig. [4] as dashed lines, are reported in Table[4]in the form 
log T e — a + b log(L/ L ), where cr is the least square estimate of 
the intrinsic dispersion. Figures [4] (a)-(f) show the evolutionary 
tracks superimposed with the quadratic fits of the edges of the 
fundamental mode instability strips and the linear fits of the first 
overtone ones for the cases I,II,III,IV,V,VI of TableQ] 

To analyse the effect of the chemical composition uncer- 
tainty, the upper panel of Fig.[3]shows the instability strips com- 
puted for the different metallicities corresponding to the cases I, 
II, III of Table [TJ whereas in the lower panel of the same figure 
the effect on the instability strip of the Y variation within its un- 
certainty is shown (case s I and IV of Table [TJl- As already found 
in our previous papers dBono et al.lll999bt ICaputo et al.l feOOO) 
the predicted instability strip becomes redder as the metallicity 
increases, while increasing the helium content from 0.26 to 0.28 
at Z = 0.005 only produces a slight narrowing of the instability 
strip at the highest luminosity levels, likely due to the slightly re- 
duced efficiency of the H ionization region in driving pulsation. 

5.3. The mass-luminosity relation 

Inspection of Fig. [4] suggests that the loop extension generally 
increases with increasing mass, but, in some cases, it de c reases 
for the 5 M model. As discussed by Castell an! et al.l {l990) 
this is because for masses around this value the first dredge-up 
(which is active for lower masses) becomes inefficient and the 
star reaches the central He burning phase with a reduced exter- 
nal He abundance (and thus with a reduced loop extension) with 
respect to the lower masses. 

For each model set we found the minimum (M m ,„) and max- 
imum (M„ Mt ) mass which cross the instability strip. The M mflt 
and M„„„ values are listed in Table [5] except for the cases in- 
cluding overshooting. The minimum mass in general increases 
with the metallicity and the helium content while the maximum 
mass increases for a decrease of the metallicity and an increase 
of Y. 

To derive the ML relation (see e.g. Bono et al. 2000) we take 
into account the three crossings of the instability strip: the first 
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Case 


a 






b 




c 




a 


Fundamental blue edge 


I (std) 


3.744 


± 


0.020 


0.039 


± 0.012 


-0.009 


± 0.002 


0.002 


II (Z/„w) 


3.751 


+ 


0.058 


0.028 


± 0.032 


-0.006 


± 0.004 


0.005 


III (Z high ) 


3.594 


+ 


0.060 


0.137 


± 0.037 


-0.024 


± 0.005 


0.005 


IV (Y high ) 


3.672 


+ 


0.032 


0.086 


± 0.018 


-0.016 


± 0.003 


0.003 


V (Yhigln Zhigh) 


3.748 


+ 


0.104 


0.046 


± 0.063 


-0.011 


± 0.009 


0.008 


VI (Z high , GS98) 


3.599 


+ 


0.076 


0.131 


± 0.046 


-0.023 


± 0.007 


0.007 


Fundamental red edge 


I (std) 


3.878 


+ 


0.112 


-0.028 


± 0.064 


-0.006 


± 0.009 


0.010 


II (Z to ) 


4.071 


+ 


0.261 


-0.146 


± 0.147 


0.012 


± 0.020 


0.025 


HI (Z Mgh ) 


3.728 


+ 


0.059 


0.071 


± 0.036 


-0.023 


± 0.005 


0.004 


IV (Y high ) 


3.988 


+ 


0.118 


-0.098 


± 0.068 


0.004 


± 0.009 


0.010 


V (Y high , Z high ) 


3.788 


+ 


0.032 


0.031 


+ 0.019 


-0.016 


± 0.003 


0.002 


VI (Z Ugh , GS98) 


3.731 


+ 


0.031 


0.064 


± 0.019 


-0.020 


± 0.003 


0.003 




4.4 4.2 4 3.8 3.6 4.4 4.2 4 3.S 3.6 

log T t Log T. 



Fig. 4. Panels (a) to (f): Evolutionary tracks, quadratic fits of the fundamental mode instability strip boundaries and linear fits of the 
first overtone ones, for the cases I to VI of Table [TJ 
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FBE quadratic fit 



FRE quadratic fit 



Table 8. Maximum (P„ raA ) minimum (P mi „) periods in days de- 
rived from the synthetic Cepheid distributions. 




— z = 


0.005 


— • z= 


0.008 


— z = 


0.0035 



3.65 
log T e 



FBE quadratic fit 



FRE quadratic fit 



3«: 




365 
l°gT e 



Fig. 3. Upper panel: effect on the fundamental instability strip 
edges of a metallicity variation within its uncertainty (cases 
I,H,III of Table [TJ. Lower panel: as in the upper panel but 
for a variation of the original helium abundance (cases I, IV 
of Table [TJ. In both panels the quadratic fits corresponding to 
Table[3]are shown. 

Table 5. Maximum (M max ) and minimum (M m ,„) mass entering 
the instability strip, for the various model sets. 



Case 


M mi „ (Ma) 


M„ MA . (Mo) 


I (std) 


4 


12 


n (Z low ) 


3 


14 


HI (Z high ) 


4 


8 


IV (Y high ) 


4 


13 



crossing during the H burning and the second and third cross- 
ings during the He burning phases. The mean luminosity is ob- 
tained by averaging the luminosity behaviour during the three 
subsequent crossings with the corresponding evolutionary times. 
Following this procedure, for each case from I to IV of Table [TJ 
we derived the corresponding ML relation (the results are avail- 
able on line in Table [6). 

Since one of the purposes of this work is to investigate the 
spread of the ML relation due to uncertainties on the chemical 
composition, we calculated the global ML relation by fitting to- 
gether all the ML data sets (cases I-IV of Table [TJ. As a result 
we find: 

log(L/L ) = 1.015(+0.021) + 3.150(+0.024) log(M/M ) (3) 

with a standard deviation of cr — 0.025; the errors include the 
effect of varying the chemical abundances among the different 
cases. The comparison with the ML relations obtained for the 
individual cases separately (see Table [6] on line) suggest that 



Case 


P m ,„ (days) 


P„„„ (days) 


I (std) 


2.5 


72.9 


II (ZfoH,) 


1.4 


87.6 


III (Zhigh) 


2.1 


21.9 


IV (Yhigi,) 


2.6 


92.0 



the coefficients of the global ML relation are consistent with 
the ones of the individual relations apart from case III (Z;,,-^) 
that corresponds to a significantly steeper relation, with a fainter 
zero-point, than all the others. This occurrence produces an in- 
crease of the intrinsic dispersion when passing from the individ- 
ual ML relations to the global one. 

To compare the present result with a previous published rela- 
tion based on evolutionary and pulsation models including older 
physical inputs, we com pared the presen t ML relation, given by 
Eq.[3] with the one bv lBono et al.l (l2000b . which explicitly takes 
into account the dependence on the chemical composition, com- 
puted for the values of Y and Z corresponding to our standard 
case: T=0.26, Z=0.005. We found tha t the M L relation derived 
in this paper is flatter than lBono et al.l (2000) with the disagree- 
ment reaching about 0.1 dex at masses around 10 M . However 
for M< ~ 8M Q the two relations are consistent within the uncer- 
tainties. 



5.4. Period-Luminosity relations 

For the cases I to IV of TableQTwe calculated the PLT e relations 
for the fundamental and first overtone models; the tables with the 
coefficients for these relations are available on line in Table [71 

We derived the synthetic PL relations by populating the 
predicted instability strip according to an assumed mass distri- 
bution and adopting, for each mass and chemical composition, 
the corresponding ML relation derived on the basis of the 
corresponding set of evolutionary tracks (see Subsec lB3b . More 
in detail, for each treated case (I to IV of Table [TJ, we built 
a synthetic population by extracting masses according to a 
Salpeter initial mass function (IMF) (oc M~ 235 ). The range 
of generated masses is given by the values of M mfl v and M„„„ 
for the cases I to IV listed in Table [5] In agreement with the 
procedure adopted in previous papers dCaputo et al.l [2000; 
iKennicuttet al.ll 19981) . we fixed to JV=1000 the total number of 
extractions falling into the instability strip. To each synthetic 
star we assigned the luminosity derived from the corresponding 
ML relation. Once the luminosity is fixed, from the quadratic fits 
of the fundamental edges, derived in Subsec |BT2l we obtained 
the corresponding temperature of the FRE and the FBE. In the 
temperature range from the T e of the FBE (decreased by 50 K) 
to the T e of the FRE (increased by 50 K), we then extracted 
randomly the temperature T e of the synthetic star. For each 
extracted T e (in K) the period (in days) is obtained from the 
P = f(L, M, T e ) relation, resulting from a linear regression 
through the fundamental models of each sel ected case. Then, b y 
means of the static model atmospheres by Bess el et all 0998), 
for the explored chemical compositions, we converted the 
model intrinsic luminosities into the various photometric bands, 
obtaining synthetic distributions in period- magnitude planes. 
The PL relations are derived by means of linear (or quadratic) 
regressions through these model distributions. 
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Table 9. Theoretical linear period-luminosity relations in differ- 
ent photometric bands for fundamental pulsators in the standard 
and global cases. 



Case 


a 


b 




O" 




std 


-1.04 ±0.03 


-2.30 ± 


0.03 


0.32 


all 


-0.87 ± 0.01 


-2.43 ± 


0.01 


0.32 


std cut 


-0.95 ± 0.03 


-2.42 ± 


0.03 


0.30 


all cut 


-0.83 ± 0.01 


-2.49 ± 


0.02 


0.30 


M v 


std 


-1.40 ± 0.02 


-2.69 ± 


0.02 


0.23 


all 


-1.28 ±0.01 


-2.77 ± 


0.01 


0.23 


std cut 


-1.34 ± 0.02 


-2.77 ± 0.02 


0.22 


all cut 


-1.25 ± 0.01 


-2.82 ± 


0.01 


0.22 


M R 


std 


-1.65 ±0.02 


-2.83 ± 


0.02 


0.20 


all 


-1.55 ±0.01 


-2.91 ± 


0.01 


0.20 


std cut 


-1.59 ±0.02 


-2.91 ± 


0.02 


0.18 


all cut 


-1.52 ±0.01 


-2.95 ± 


0.01 


0.19 


Mi 


std 


-1.92 ±0.01 


-2.96 ± 


0.02 


0.17 


all 


-1.84 ± 0.01 


-3.01 ± 


0.01 


0.17 


std cut 


-1.87 ±0.02 


-3.02 ± 0.02 


0.15 


all cut 


-1.82 ±0.01 


-3.05 ± 


0.01 


0.16 


Mj 


std 


-2.238 ± 0.010 


-3.149 ± 0.011 


0.12 


all 


-2.184 ± 0.005 


-3.184: 


± 0.005 


0.12 


std cut 


-2.206 ±0.011 


-3.191 


± 0.012 


0.11 


all cut 


-2.166 ± 0.005 


-3.210 : 


± 0.006 


0.11 


M K 


std 


-2.499 ± 0.006 


-3.329 : 


± 0.006 


0.07 


all 


-2.470 ± 0.003 


-3.343 : 


± 0.003 


0.08 


std cut 


-2.481 ± 0.007 


-3.353 : 


± 0.008 


0.07 


all cut 


-2.459 ± 0.003 


-3.359 : 


± 0.004 


0.07 



The maximum and minimum periods corresponding to 
the obtained synthetic distributions are reported in Table [8] In 
Tables l9l and [TTI we report the linear (Mi = a + b log P) and 
quadratic (Ma - a + b log P + c (log P) 2 ) period-luminosity 
relations for fundamental pulsators in the different photometric 
bands. In these Tables we listed the fit coefficients for case 
I of Table Q] (std) and the global fit of cases I to IV taken 
together (all), as well as the results for the PL relations obtained 
by rejecting periods longer than logP > 1.5 (cut). The case 
all provides an estimate of the effect of the uncertainties on 
the chemical abundances on the predicted PL relations. The 
linear and quadratic PLx relations, computed for the single cases 
II, III and IV of TableQ] are available on line in Tables [T0land[T2l 

Passing from the individual cases to the global PL re- 
lations, the coefficients and the intrinsic dispersion do not 
vary significantly because the effects of these variations in the 
abundances are smaller than the effect on the intrinsic scatter 
due to the finite width of the instability strip, especially in the 
optical filters. Moreover, inspection of Table [9] suggests that 
the PL relations become more linear, less dependent on the 
chemical composition and with a smaller intrinsic dispersion 
passing from the optical to the near-infrared filters, c onfirming 
previous observational ( Madore & Freedman 1991) and the- 
oretical dBono et all Il999bt ICaputo et al.l |2000|) results. The 
PL relations obtained for log P < 1.5 are derived because, as 



shown in lCaputo et al.l d2000l) . removing the longest periods, the 
period-magnitude distribution is better reproduced by a linear 
relation also in the optical bands. Moreover, many observed 
Cepheid samples cover this period range. 

The deviation f rom linearity of C e pheid PL relations 
is a debated issue (fSandage etall 12004 iKanbur et ail |2004 
Ngeo w et al.l 120051 : IKanbur et alJl2006l) and for t his reason we 
also co mputed quadratic PL relations. In particular. lNgeow et al.l 
d2005l) performed a statistical investigation of the LMC classi- 
cal Cepheid sample obtained from the MACHO database and 
found that the observed behaviour in period-magnitude diagrams 
is best reproduced by two linear relations, with a break at 10 
days. Nonlinear pulsation models (see also ICaputo et alJ [2000: 
iFiorentino et al.|[2002k iMarconi et al.l l2005) suggest a quadratic 
form of PL relations, in particular in the optical bands (see Table 
rm >. However the effect of this nonlinearity on the calibration 
of the extragalactic distance scale is very small (see iKoen et ail 
120071) . In Fig. [5] we plotted the predicted synthetic PL relations 
in the V (upper panel) and I (lower panel) Johnson-Cousins pho- 
tometric bands with the quadratic fits superimposed (solid lines) 
for cases I to IV taken together (case all in Table [TTI). As a com- 
parison the two linear regressions (for log P < 1 .0 and log P > 

1.0, d ashed lines) obtained observationally by ISandage et al.l 
d2004) are also shown. In panel a of Fig.|6]we show the residuals 
of the present linear PL relations for log P < 1.5 (case all cut of 
Tabled in the V, I bands and over the whole peri od range for the 
K ban d with respect to th e empirical relations b y lFreedman et al.l 
d2001l) (for V and I) and lPersson etafl d2004l) (for the K band). 
According to this plot, we find differences smaller than 0.2 mag 
(in absolute value) in all the bands. Panel b of Fig. [6] shows the 
differences between the present linear solutions with log P < 1.5 
and the theoretical PL linear relations with logP < 1.5 by 
Caput o et al.1 d2000l) . A gain we find differences within 0.2 mag 
with an opposite trend as a function of the period with respect 
to the previous comparison, thus suggesting that the the magni- 
tudes predicted by the new PL relations obtained in this paper 
are somewhat intermediate between the ones based on emp irical 
relations and the theoretical relations bv lCaputo et alj d2000h . 

As far the physical assumptions of our pulsation models are 
concerned, we have already shown that changes in the adopted 
equati on of state do not affect significantly the pulsation sce- 
nario dPetroni| [2003). On the other hand, variations of the mix- 
ing lengt h parameters used to close t he nonlinear equation sys- 
tem (see Bono & Stelling werfl 1 1 994) have been found to pro- 
duce negligible effect s on the predicted Cepheid relations (see 
Fiorentino et al .1 120071 for details). 

6. Noncanonical models 

6.1. Overshooting 

The extension of the mixing into a region that is stable fol- 
lowing the Schwar zschild criterion is generally called "over- 
shooting" (see e.g. Co rdier et al.l l2002. and references therein). 
An additional extension of the canonical convective regions 
could also be due to r otationally induced mixing (see e.g. 
iMevnet & Maeden uOOO): both phenomena can be modeled, at 
least as far as the extension of the mixing zones is concerned, 
with the same "overshooting" formalism. A complete under- 
standing of these extra-mixing mechanisms is still lacking in 
the literature due to the complexity of the involved phenomena, 
that also prevents an analyti cal approach (see e.g. IStraka et~ai1 
120051: iDem arque et al. 2007). Three-dimensional hydrodynam- 
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Table 11. Theoretical quadratic period-luminosity relations for 
fundamental pulsators in the standard and global cases. 



Case 


a 




b 


c 


cr 


M B 


std 


-0.43 


± 0.08 


-3.69 ±0.17 


0.69 ± 0.08 


0.31 


all 


-0.63 


±0.02 


-3.07 ± 0.06 


0.35 ± 0.03 


0.32 


My 


std 


-0.98 


± 0.06 


-3.66 ±0.12 


0.48 ± 0.06 


0.22 


all 


-1.10 


± 0.02 


-3.27 ± 0.04 


0.27 ± 0.02 


0.23 


M R 


std 


-1.29 


± 0.05 


-3.65 ±0.10 


0.41 ± 0.05 


0.19 


all 


-1.39 


± 0.02 


-3.33 ± 0.04 


0.23 ± 0.02 


0.19 


M, 


std 


-1.61 


±0.04 


-3.66 ± 0.09 


0.35 ± 0.04 


0.16 


all 


-1.70 


±0.01 


-3.38 ±0.03 


0.20 ± 0.02 


0.16 


Mj 


std 


-2.03 


± 0.03 


-3.63 ± 0.06 


0.24 ± 0.03 


0.12 


all 


-2.08 


±0.01 


-3.45 ± 0.02 


0.14 ±0.01 


0.12 


M K 


std 


-2.38 


±0.02 


-3.61 ± 0.04 


0.14 ±0.02 


0.07 


all 


-2.41 


±0.01 


-3.50 ± 0.01 


0.09 ± 0.01 


0.07 



ical calculations are still un der development (Eg gleton et alJ 
l2007t iMeakin & Arnettll2007b . Thus a semi-empirical approach 
is currently used by defining the extension of any additional mix- 
ing region, starting from the border of the canonical convective 
zones, in terms of the pressure scale height H v : l m = (3H V . The 
overshooting efficiency is generally estimated through the com- 
parison between theory and observations for diffe rent relevant 
quantities, and is still widely debated point (see e.g.lClaretl2007t 
iBrocato et alll2003l:lBarmina et alj|2002l IChiosi et al.ll 19921 and 
references therein). 

We calculated additional sets of evolutionary models in- 
cluding overshooting during the central H burning phase (Main 
Sequence, MS) with l m = 0.1 H p and / ov = Q.25 H B , within the 
range suggested in the recent literature (see e.g. Girard fet al.l 
2000). During the central He burning phase our c ode implements 
the cl assical semi-convection criterion (see e.g. Castellan i et al.1 
119851) at the border of the Schwarzschild convective core (which 
is driven by mechanical overshooting at the convective core 
boundary), while breathing pulses (see e.g. [Castellan i et al.l 
[1981 15 orman & Roodll 1993b are suppre ssed following the pro- 
cedure described bv lCaputo et al.l d 1989b . In the present models 
the overshooting region is fully homogenized but the tempera- 
ture gradient is kept at its radiative value. 

The inclusion of overshooting during central H burning in 
standard models leads to a higher luminosity and longer lifetime 
during the MS phase, larger He cores and luminosities and re- 
duced lifetimes during the following central He burning phase. 
This reflects in the adoption of a different mass-luminosity re- 
lation for Cepheid pulsational models. The adopted luminosity 
levels for noncanonical models in this paper are taken as brighter 
th an the canonical on es by 0.25 dex , following the prescriptions 
bv lChiosi et al.l d!993l) . We checked that it is true for current evo- 
lutionary models too. 

When the overshooting phenomenon is included, the blue 
loops in the color-magnitude diagram are less extended than in 
the absence of overshooting in such a way that for some chem- 
ical compositions intermediate-mass models do not cross the 
classical Cepheid instability strip. To avoid this problem, some 




LogP 



Fig. 5. Synthetic period-luminosity relations in the V (upper 
panel) and I (lower panel) Johnson-Cousins photometric bands 
with the quadratic fits (solid lines) corresponding to the case all 
of Table [TTI Symbols: std case (case I, black dots); Z/ 0M , (case II, 
green dots); Z/„^/, (case EH, blue dots); Yy^gh (case IV, red dots). 



authors include an "ad hoc" u ndershooting at the l ower boundary 
of convective envelopes (see iGirardi et al.ll2000l) . even if there 
is no clear observationa l evidence for this mechanism (see e.g. 
iRenzini & Ritossalll994l:rRiello et alj|2003l) . 

Relevant pulsation properties, such as the morphology and 
the amplitude of light and radial velocity curves, are affected by 
the assumed ML relation dBono et al.l l2002; Nat ale et al.l l2008). 
and, in t urn, by the amount of overshoo ting accounted for in the 
models (Caput o et alj2005l:lKellej2 008. and references therein). 
Independent theoretical and observational evidence suggests that 
the Cepheid pulsation masses are lower than the canonical evo- 
lutionary values but the amount of this discrepancy is not firmly 
established even if several au t hors find a value of the order of 10- 
15 % (e.g.lBono et alj|2001t iBeaulieu et afll200U ICaputo et al.1 
2005; lNataleet al.ll2 008. and references therein). 

To quantify the effect of overshooting, we computed models 
with standard chemical composition (Z = 0.005, Y = 0.26) for 
two selected values of the overshooting efficiency (j8=0.1 and 
/J=0.25, models XI and XIII of Table Q). However, due to the 
reduced loop extension, some intermediate mass models with 
/3-0.25 do not cross the instability strip, in disagreement with 
observations of LMC Cepheids in this region. To avoid this 
problem without including overshooting for the convective en- 
velopes, taking into account the present uncertainty in the orig- 
inal helium abundance, we calculated overshooting models for 
Z = 0.005, Y = 0.28 (model XII of Table Q) which, as already 
discussed, show a more extended loop than the standard models. 
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0.5 1 

logP 



Fig. 6. a) Residuals of the present linear PL relations for log P < 
1.5 in the V, I bands and over the whole period range for the 
K band with respect to the empirical relations; b) Differences 
between the presen t linear solutions w ith logP < 1.5 and the 
theoretical ones by Caput o et"aT] ((2000). 



Table 13. Quadratic analytical relations for the fundamental 
edges, in the non canonical standard case. 



a 




b 


c 




a 


Blue Edge 


3.904 


± 0.046 


-0.034 ± 0.027 


-0.0001 


± 0.0038 


0.002 


Red Edge 


3.816 


+ 0.078 


0.008 ± 0.046 


-0.0124 


± 0.0065 


0.004 



In this way, even models with /3=0.25 populate the instability 
strip. 

In Table [T3l the coefficients of the quadratic fits of the FRE 
and FBE for models with Z ov = 0.25// p and standard chemical 
composition are reported. 

The upper panel of Fig. [7] shows the evolutionary tracks with 
Zov = 0.25H p , Z=0.005, 7=0.28 and the quadratic fit of the non 
canonical strip, while the lower panel of the same figure shows 
the evolutionary tracks with / ov = 0.1H p , Z=0.005, T=0.26 and 
the quadratic fits of the canonical (dot-dashed lines) and non 
canonical (solid lines) instability strip boundaries calculated for 
standard chemical composition. We notice that the case with 
j8=0. 1 is not significantly different from the canonical one and 
for this reason in the following we will concentrate on the case 
with £=0.25. 

The ML obtained with the set of tracks with / ov = 0.25// p , 
Z=0.005, T=0.28 combined with the noncanonical strip of 
Fig. [7] upper panel, is the following: 

log(L/L G ) = 1.260(+0.061) + 3.187(+0.076) log(M/M ) (4) 




4.4 4.2 4 3.8 3.6 

togT 



Fig. 7. Evolutionary tracks with F=0.28, Z=0.005, l m = 0.25H p 
and the quadratic fit of the noncanonical instability strip bound- 
aries (upper panel); evolutionary tracks for T=0.26, Z=0.005, 
/ ov = 0.1 H p and the quadratic fit of the canonical (dot-dashed 
lines) and noncanonical (solid lines) instability strip boundaries 
for standard chemical composition (lower panel). 



with <x = 0.022. As expected this ML relation has a brighter 
zero-point than the canonical one (see Eq.|3) by about 0.25 dex, 
and the same slope within the errors. This difference can be also 
noted from the comparison shown in Fig. [8] 

For the chemical composition of case I, in Table[TJwe calcu- 
lated the PLT e relations for noncanonical fundamental and first 
overtone models; the tables with the coefficients for these rela- 
tions are available on line in Table [T4l 

Adopting the same procedure as for canonical models we 
derived the multi-band period-luminosity relations for [3 = 0.25, 
Z=0.005, 7=0.28. The coefficients of the linear PL in the V, I, K 
bands are reported in Table 1 171 whereas the linear and quadratic 
PL relations in the other photometric bands are available on line, 
in Tables [BJ and [16] 

A comparison between the canonical (case std of Table |9]l 
and the noncanonical PL linear relation in the V, I, K bands (see 
Table [IT) is shown in Fig. [9] where the errorbars refer to the in- 
trinsic dispersions of the canonical relations. We notice that both 
the zero point and the slope of the PL relations are affected by 
the inclusion of mild overshooting. The noncanonical V band re- 
lation predicts brighter magnitudes than the canonical ones in the 
very short period range and fainter magnitudes than the canoni- 
cal ones for periods longer than about 5 days, whereas the non- 
canonical I and K band PL predict almost identical magnitudes 
to the canonical ones at the shortest periods and fainter magni- 
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Fig. 8. Comparison between the mass-luminosity relation from 
Eq. El (7=0.28, Z=0.005, l ow = 0.25# p ) and the global ML rela- 
tion from Eq.[3] 



Table 17. Theoretical PL relations for fundamental noncanoni- 
cal pulsators (see text). 



Case 


a 




b 






M v 


ncOS 
ncOS cut 


-1.60 
-1.47 


± 0.04 
± 0.04 


-2.34 
-2.47 


±0.03 
±0.04 


0.29 
0.28 


M, 


ncOS 
ncOS cut 


-2.01 
-1.92 


± 0.03 
± 0.03 


-2.7 i 
-2.8 i 


0.02 
0.03 


0.21 
0.20 


M K 


ncOS 
ncOS cut 


-2.43 
-2.39 


± 0.01 
± 0.01 


-3.21 
-3.25 


±0.01 
±0.01 


0.08 
0.08 



Table 18. M min , M max , (as in Table |5] |5]l for Z Q 

0.25// n 



Case 




M 


p 

1 mm 


p 

L max 




(M G ) 


(Mo) 


(days) 


(days) 


XII (l m = 0.25/fp) 


4 


9 


4.1 


72.4 



-8 



-4 



-8 
-6 

s" 

-4 



-8 





1 1 1 1 1 1 1 1 1 1 1 






_ _ std case, non canonical (Z=0.005, Y=0.26, l m =0.25 H ) 
std case, canonical (Z=0.005, Y=0.26) 
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Fig. 9. Comparison, in the V, I and K photometric bands, be- 
tween the present canonical and noncanonical linear PL relations 
(solid and dashed lines). 



tudes than their canonical counterparts elsewhere, with the effect 
increasing toward longer periods. At a period of 10 days, the ef- 
fect of applying the noncanonical relations instead of the canon- 
ical ones would be of about 0.15, 0.17 and 0.19 mag in the V, I 
and K band respectively, whereas at 30 days it would be as large 
as about 0.32, 0.29 and 0.24 in the V, I and K band respectively. 

Referring to the upper panel of Fig. [7] we report in Table [T8l 
the minimum and maximum masses with blue loops that cross 
the red edge of the instability strip, together with the correspond- 
ing values of the minimum and maximum period. 

6.2. Mass loss 

The possibility that mass loss might solve the long-standing 
problem of the mass discrepancy of Cepheids between the 
value inferred from pulsational and evoluti onary models has 
been suggested by various authors (see e.g. iBono et ail 120021; 
ICaputo et al J 12005). The argument is that in order to solve the 



problem - that is, to provide higher luminosity for He-burning 
stars of a given mass - a larger He-core mass is needed. As pre- 
viously discussed, this is precisely the effect of the convective 
core-overshooting. On the other hand, a similar result can be 
mimicked by a standard model affected by significant mass loss. 
In order to check this, we computed the evolution of a 5 M Q 
star taking into account the effect of a mild core overshooting 
(Z ov = 0.25//p) during the H-burning phase. We attempted to re- 
produce the minimum luminosity (LS ' ) of this model during the 
central He-burning by standard models affected by mass loss. 
The mass loss process is not yet fully understood and it is still 
lacking satisfactory knowledge of its efficiency along the various 
evolutionary phases. As a consequence, we are forced to adopt 
different mass loss prescriptions in the following numerical ex- 
periments in order to study the problem. Thus we computed the 
evolution of stars with mass in the range between 6.0 and 7 
M Q with different mass loss prescriptions adopting the classical 
Schwarzschild criterion to determine the convective core bound- 
ary; the results, discussed below, are shown in Fig. [10] The mod- 
els were evolved with constant mass during the main sequence 
and subgiant phase, while a substantial mass loss, with a con- 
stant rate, was imposed during the red giant phase. Once the to- 
tal mass of the models reached the desired value, the subsequent 
evolution was performed at constant mass. Following such a pro- 
cedure, we computed for example the evolution of the central 
He-burning phase of a 5.3 M Q with a progenitor star of 6.2 M , 
whose behavior in the HR-diagram is like that of our target, the 5 
M Q with core overshooting. However, further numerical experi- 
ments showed that the result is not unique; models with different 
initial masses and different prescriptions for the mass loss have 
a minimum luminosity during the central He-burning very sim- 
ilar to . Such a result was expected since the Cepheid ML 
relation depends appreciably on both the mass-loss rate and the 
evolutionary phase at which it is turned on. An example is given 
by a star with initial mass of 6.1 M evolved from the ZAMS, 
with the classical Reimers mass loss rate (with 77=1); even if 
this model describes a much more extended loop in color in the 
HR-diagram, the minimum luminosity of the central He burning 
phase is the same. The target luminosity can also be reproduced 
by a more extreme model: a 3 M with a progenitor of 6.4 M . In 
this case, the color of the He-burning model is significantly red- 
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logT 

Fig. 10. Comparison of the minimum luminosity during the cen- 
tral He burning phase for a 5 M Q model with l m = Q.25H p and 
models of different original masses with different massloss pre- 
scriptions. 

der than the model with core overshooting. These few numerical 
examples prove that it does not make sense to try to compute a 
ML relationship for standard models with mass loss, unless a pri- 
ori mass loss prescriptions have been fixed. Moreover, the mass- 
loss solution of the Cepheid mass discrepancy sounds rather ad 
hoc, at least until a physical process responsible for the huge 
required mass loss is found and understood. 

7. Conclusion 

On the basis of an updated set of evolutionary and pulsation 
models for Classical Cepheids in the LMC, we have investigated 
the effect of the uncertainties on the chemical composition and 
on the physical assumptions adopted in the codes on the most 
relevant pulsation observables and thus on the theoretical cal- 
ibration of the Cepheid distance scale. We found that present 
uncertainties on relevant nuclear reaction rates have only a neg- 
ligible effect on evolutionary and pulsational theoretical predic- 
tion. Moreover the uncertainties in the metal and helium abun- 
dances affect the results of the evolutionary computations but 
do not significantly changethe pulsation scenario. On the other 
hand, the still present uncertainties on the efficiency of the over- 
shooting phenomenon in the previous H burning phase and on 
the mass-loss rates are found to be the most important source 
of uncertainty in the theoretical Cepheid mass-luminosity and 
period-luminosity relations. In particular, using a theoretical PL 
relation that relies on the assumption of mildly overshooting 
evolutionary models, one infers distances that are significantly 
shorter than the values obtained when a canonical theoretical PL 
is used, especially in the optical bands. Therefore we conclude 
that the uncertainty on the Cepheid ML relation is expected to 
significantly affect the Cepheid calibration of the extragalactic 
distance scale and in turn the evaluation of the Hubble constant. 
The application of the theoretical analysis performed in this pa- 
per to LMC Cepheid data will be addressed in a forthcoming 
paper (Valle et al. 2009 in preparation). 
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Fig. 1. Evolutionary tracks with standard chemical composition 
with three different values of the 12 C(q','>') 16 rate (see text). 



Table 2. Coefficients of the linear analytical relations for the fun- 
damental red and blue edges. 



Case 


a 




b 






a 


Blue Edge 


I (std) 


3.848 


± 0.010 


-0.022 


± 


0.002 


0.005 


II (Z low ) 


3.828 


± 0.013 


-0.016 


± 


0.003 


0.006 


HI (Z«g*) 


3.853 


± 0.029 


-0.025 


± 


0.008 


0.012 


IV {Y Mgh ) 


3.866 


± 0.020 


-0.026 


± 


0.005 


0.009 


V (Yingh, Z/„- g /,) 


3.870 


± 0.020 


-0.029 


± 


0.006 


0.008 


VI (Z kigk , GS98) 


3.858 


± 0.022 


-0.027 


± 


0.007 


0.011 


Red Edge 


I {std) 


3.955 


± 0.020 


-0.073 


± 


0.005 


0.010 


II (Zw) 


3.922 


± 0.048 


-0.061 


± 


0.012 


0.023 


in (Z kigk ) 


3.974 


± 0.028 


-0.082 


± 


0.008 


0.012 


IV (Y high ) 


3.934 


± 0.021 


-0.066 


± 


0.005 


0.009 


V (Yhigh, Zhigh) 


3.965 


± 0.019 


-0.078 


± 


0.005 


0.008 


VI (Z Ugh , GS98) 


3.958 


± 0.018 


-0.074 


± 


0.005 


0.009 



Table 6. ML relations for the selected cases; log(L/L ) = a + 
b log(M/M ). 



Case 


a 




b 




a 


I (std) 


1.000 


± 0.017 


3.160 


± 0.019 


0.009 


II CZtoJ 


1.040 


± 0.010 


3.110 


± 0.011 


0.008 


HI (Z Ugh ) 


0.736 


± 0.008 


3.455 


± 0.010 


0.002 


IV (Y high ) 


1.095 


± 0.043 


3.074 


± 0.047 


0.024 



Table 10. Theoretical linear period-luminosity relations in dif- 
ferent photometric bands for individual cases. 



Case 


a 


b 


<x 


M B 


Zlow 


-0.72 ± 0.02 


-2.61 ± 0.02 


0.27 


Zhigh 


-1.00 ± 0.03 


-2.18 ± 0.04 


0.29 


yhigh 


-1.03 ± 0.03 


-2.28 + 0.03 


0.37 


My 


Zlow 


-1.16 ± 0.01 


-2.90 ± 0.01 


0.20 


Zhigh 


-1.40 ± 0.02 


-2.61 ± 0.03 


0.21 


Yhigh 


-1.39 ± 0.02 


-2.68 ± 0.02 


0.27 


M R 


Zlow 


-1.44 ± 0.01 


-3.00 ± 0.01 


0.16 


Zhigh 


-1.66 ± 0.02 


-2.77 ± 0.02 


0.18 


Yhigh 


-1.63 + 0.02 


-2.82 + 0.02 


0.23 


M, 


Zlow 


-1.75 ± 0.01 


-3.10 + 0.01 


0.14 


Zhigh 


-1.93 ± 0.02 


-2.90 ± 0.02 


0.15 


Y/iigh 


-1.90 ± 0.02 


-2.95 ± 0.02 


0.20 


Mj 


Zlow 


-2.117 ± 0.006 


-3.240 + 0.007 


0.10 


Zhigh 


-2.247 ± 0.011 


-3.119 ± 0.014 


0.11 


Yhigh 


-2.217 + 0.013 


-3.140 ±0.013 


0.14 


M K 


Zlow 


-2.430 + 0.003 


-3.373 ± 0.004 


0.06 


Zhigh 


-2.507 ± 0.007 


-3.317 + 0.009 


0.07 


Yhigh 


-2.477 ± 0.008 


-3.320 ± 0.008 


0.09 



Table 12. Theoretical quadratic period-luminosity relations for 
fundamental pulsators in the various individual cases. 



Case 


a 


b 




c 








Z{ow 
Zhigh 
Yhigh 


-0.62 + 0.03 
-0.25 ± 0.08 
-0.74 ± 0.10 


-2.92 
-4.32 
-2.92 


± 0.08 
± 0.22 
± 0.20 


0.18 
1.35 
0.30 


± 0.04 
± 0.14 
± 0.09 


0.27 
0.27 
0.36 


My 


Zlow 
Zhigh 
Yhigh 


-1.10 ± 0.02 
-0.87 ± 0.06 
-1.19 ± 0.07 


-3.08 
-4.12 
-3.09 


± 0.06 
± 0.16 
± 0.14 


0.11 
0.96 
0.20 


± 0.03 
± 0.10 
± 0.07 


0.20 
0.20 
0.27 



Mr 

Zi„w -1.40 + 0.02 -3.15 ± 0.05 0.08 ± 0.03 0.16 

Zhigh -1.21 ± 0.05 -4.05 ±0.14 0.81 ± 0.09 0.17 

Yhigh -1.47 ±0.06 -3.17 ±0.12 0.16 ± 0.06 0.23 

~M, 

Ziow -1.71 ± 0.02 -3.22 ± 0.04 0.07 ± 0.02 0.14 

Z Hg h -1.54 ±0.04 -4.01 ± 0.12 0.70 ± 0.07 0.14 

Yhigh -1.76 ± 0.05 -3.25 ±0.10 0.14 ± 0.05 0.19 

M7 

Ziow -2.09 ±0.01 -3.32 ± 0.03 0.04 ± 0.02 0.10 

Z Hgh -1.97 ± 0.03 -3.90 ± 0.08 0.50 ± 0.05 0.10 

Yhigh -2.13 ± 0.04 -3.34 ± 0.08 0.09 ± 0.04 0.14 

Ziow -2.42 ± 0.01 -3.41 ± 0.02 0.02 ± 0.01 0.06 

Z U gh -2.34 ± 0.02 -3.80 ± 0.05 0.31 ± 0.03 0.07 

Yhigh -2.43 ± 0.02 -3.43 ± 0.05 0.05 ± 0.02 0.09 
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Table 7. Canonical PLT e relations for the selected cases; logP = a + b log(L/L ) + c (log T e ). 



Case 


a 




b 




c 




<X 


Fundamental - CANONICAL 


I (std) 


10.725 


± 0.180 


0.687 


± 0.004 


-3.261 


± 0.046 


0.013 


II (Zw) 


10.475 


± 0.320 


0.696 


± 0.005 


-3.202 


± 0.082 


0.016 


in (Z high ) 


10.908 


+ 0.173 


0.677 


± 0.004 


-3.303 


± 0.044 


0.009 


IV (Y high ) 


11.265 


± 0.270 


0.681 


± 0.005 


-3.398 


± 0.069 


0.013 


First Overtone - CANONICAL 


I (std) 


10.827 


± 0.393 


0.633 


± 0.004 


-3.282 


± 0.101 


0.003 


II (Z w ) 


10.826 


±0.187 


0.631 


± 0.002 


-3.279 


± 0.048 


0.002 


in (Z high ) 


11.033 


± 0.404 


0.630 


± 0.003 


-3.337 


± 0.105 


0.004 


IV (F fcsA ) 


10.897 


± 0.347 


0.633 


± 0.003 


-3.299 


± 0.090 


0.002 



Table 14. Non canonical PLT e relations: logP = a + b log(L/L G ) + c (log T e ). 



Case 


a 




b 




c 




cr 


Fundamental - Non canonical 


I (std) 


10.290 


± 0.290 


0.701 


+ 0.005 


-3.144 


± 0.075 


0.018 


First Overtone - Non canonical 


I (std) 


10.840 


± 0.130 


0.639 


± 0.002 


-3.278 


± 0.032 


0.001 



Table 15. Theoretical linear period-luminosity relations for the 
fundamental noncanonical pulsators. 



Case 


a 


b 




M B 


ncOS 


-1.40 ± 0.05 


-1.81 ± 0.04 


0.41 


ncOS cut 


-1.22 ± 0.06 


-2.00 ± 0.06 


0.40 


M R 


ncOS 


-1.79 ±0.03 


-2.54 ± 0.03 


0.25 


ncOS cut 


-1.68 ±0.04 


-2.65 ± 0.03 


0.24 


Mj 


ncOS 


-2.24 ± 0.02 


-2.96 ± 0.02 


0.14 


ncOS cut 


-2.18 ±0.02 


-3.03 ± 0.02 


0.14 



Table 16. Theoretical quadratic period-luminosity relations for 
fundamental noncanonical pulsators. 



Case 


a 


b 


c 


cr 


M B 


ncOS 


-0.87 ± 0.20 


-2.78 ± 0.36 


0.41 ± 0.15 


0.41 


ncOS cut 


-2.45 ± 0.29 


0.42 ± 0.56 


-1.13 ±0.26 


0.39 


M v 


ncOS 


-1.23 ± 0.15 


-3.01 ± 0.26 


0.29 ±0.11 


0.29 


ncOS cut 


-2.39 ± 0.21 


-0.67 ± 0.4 


-0.84 ±0.19 


0.28 


M R 


ncOS 


-1.47 ± 0.12 


-3.12 ± 0.22 


0.25 ± 0.09 


0.25 


ncOS cut 


-2.44 ± 0.17 


-1.14 ± 0.33 


-0.7 ±0.16 


0.24 


Mi 


ncOS 


-1.72 ±0.10 


-3.22 ± 0.18 


0.22 ± 0.08 


0.21 


ncOS cut 


-2.54 ± 0.15 


-1.56 ± 0.28 


-0.58 ±0.13 


0.2 


Mj 


ncOS 


-2.07 ± 0.07 


-3.28 ± 0.13 


0.14 ± 0.05 


0.14 


ncOS cut 


-2.64 ± 0.10 


-2.13 ± 0.20 


-0.42 ± 0.09 


0.14 


M K 


ncOS 


-2.35 ± 0.04 


-3.36 ± 0.07 


0.06 ± 0.03 


0.08 


ncOS cut 


-2.66 ± 0.06 


-2.72 ±0.11 


-0.25 ± 0.05 


0.08 



